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ABSTRACT 

We examine the X-ray luminosity scaling relations of 31 nearby galaxy clusters from the Representative XMM-Newton Cluster 
Structure Survey (REXCESS). The objects are selected only in X-ray luminosity, optimally sampling the cluster luminosity function. 
Temperatures range from 2 to 9 keV, and there is no bias toward any particular morphological type. To reduce measurement scatter 
we extract pertinent values in an aperture corresponding to Rsoo, estimated using the tight correlation between Yx (the product of gas 
mass and temperature) and total mass. The data exhibit power law relations between bolometric X-ray luminosity and temperature, 
Yx and total mass, all with slopes that are significantly steeper than self-similar expectations. We examine the possible causes for the 
steepening, finding that structural variations have little eff'ect and that the primary driver appears to be a systematic variation of the gas 
content with mass. Scatter about the relations is dominated in all cases by the presence of cool cores. The natural logarithmic scatter 
about the raw X-ray luminosity-temperature relation is about 70 per cent, and about the X-ray luminosity- 7^ relation it is 40 per 
cent. Systems with more morphological substructure show similar scatter about scaling relations than clusters with less substructure, 
due to the preponderance of cool core systems in the regular cluster subsample. Cool core and morphologically disturbed systems 
occupy distinct regions in the residual space with respect to the best fitting mean relation, the former lying systematically at the high 
luminosity side, the latter lying systematically at the low luminosity side. Simple exclusion of the central regions serves to reduce the 
scatter about the scaling relations by more than a factor of two. The scatter reduces by a similar amount with the use of the central gas 
density as a third parameter. Using Yx as a total mass proxy, we derive a Malmquist bias-corrected local luminosity-mass relation and 
compare with other recent determinations. Our results indicate that luminosity can be a reliable mass proxy with controllable scatter, 
which has important implications for upcoming all-sky cluster surveys, such as those to be undertaken with Planck and eROSITA, and 
ultimately for the use of the cluster population for cosmological purposes. 

Key words. Cosmology: observations. Cosmology: dark matter. Galaxies: cluster: general, (Galaxies) Intergalactic medium. X-rays: 
galaxies: clusters 

1, Introduction allowed further investigation with increasingly better quality 

The X-ray luminosity is an observationally attractive quantity 

because of the relative ease with which it can be measured, and The density squared {nl) dependence of the X-ray emission 
thus it is a key parameter for cosmological applications of the means that luminosity measurements are very sensitive to the 
galaxy cluster population. For a fully virialised cluster formed exact physics of the gas near the cluster core. Mechanisms such 
through pure gravitational collapse, the X-ray luminosity L is as rapid radiative cooling or merging can change the thermody- 
determined solely by the mass and distribution of gas in the in- namic state of this core gas, introducing scatter into the various 
tracluster medium (ICM), and the X-ray temperature T is deter- luminosity scaling relations. Since our knowledge of the abso- 
mined by the depth of the potential well in which the ICM rests, lute extent of the scatter limits the constraints that can be put on 
Correlations between these two basic quantities were found cosmological models with the cluster population, the magnitude 
in the very early days of X-ray observations of clusters, even of the scatter, its source(s), and how to correct for it constitute 
while the thermal nature of the emission was still under debate some of the most importan t open issues in the st udy of clusters 
("Mitchell etaLl 119771: iMushotzkv et"an 119781: iHenrv & Tucked (see e.g. iLima & Hul [2QQ5h. iFabian et al.1 11994') were the first 
1979). Initial results from these works suggested that the slope to note that the oflfset of a cluster from the mean relation was 
of the luminosity temperature relation was steeper than expected connected to the presence of a cool cor e, motivating examina- 
from gravitational collapse alone. The launch of the EXOSAT tion of methods to correct for this eflTect. i MarkevitchI ri998) de- 
and Einstein observatories ena bled the first systematic studies rived quantities corrected for the presence of cooling cores by 
of large samples of clusters (Edge & Stewartl ll99ll iDavid et all exclusion of the emission from the central region and the in- 
[1993). and the subsequent launch of ROSAT, ASCA and Ginga troduction of a second spectral com ponent in the temperature 

estimation; lArnaud & EvrardI (Il999l) determined the luminosity 
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weak or non-existent cool cores. More recent efforts have aimed 
at reducing scatter by using the peak surface brightness as a third 
parameter (O'Haraetal. 2006). At the same time it had long 
been suspected that merging events also contributed to scatter 
about the mean relation. This effect has been investigated with 
increasingly sophisticated numerical simulations, which have 
shown that while major mergers can indeed boost both lumi- 
nosity and temperature, the boosting appears to be short lived 
and the net movement in the luminosity te mperature plane is 
approximately parallel to the mean relation ([Ricker & SarazinI 
l2QQll:lRitchie & Thomasll2002l: iHartlev et"aDl2008h . 

In the present paper we re-investigate the luminosity scal- 
ing relations with REXCESS (Bohringer et al. 2007), a sample 
of 33 local (z < 0.2) cluste rs drawn from the REFLEX cata- 
logue teohringer et al.ll2004l) . all of which have been observed 
with a single satellite, XMM-Newton, with the aim of minimis- 
ing effects due to instrumental cross calibration uncertainties. 
The unique sample selection strategy, in which clusters have 
been selected by luminosity only, in such a way as to have close 
to homogeneous coverage of luminosity space, delivers an op- 
timal sampling of the luminosity function of the cluster popu- 
lation with no bias towards any particular morphological type. 
Moreover, distances were optimised so that the angular scale of 
the objects is such that ~ R500 falls well within the XMM-Newton 
field of view, allowing detailed local background modelling to be 
undertaken, increasing the precision of measurements at large 
radii as compared to more nearby clusters which often fill the 
field of view. Since the basic selection criterion is X-ray lumi- 
nosity, REXCESS should be representative of any local, unbi- 
ased high-quality X-ray survey, of the type applicable to testing 
of cosmological models. 

In the following, we first use the representative nature of the 
REXCESS sample to investigate the raw luminosity scaling rela- 
tions, finding that the slopes are steeper than expected if the gas 
is heated purely by gravitational processes and that cooling cores 
are the dominant contributor to scatter about them. Dividing the 
data into subsamples, we investigate the effect of cool cores and 
morphological disturbance on a cluster's position with respect 
to the mean relation, finding that the former lie systematically at 
the high luminosity side, and the latter lie systematically at the 
low luminosity side. We then investigate two different methods 
to minimise scatter: simple exclusion of the central region, and 
use of the central gas density as a third parameter in scaling law 
fitting, finding that both methods result in a significant reduc- 
tion in the dispersion about the best fitting relations. Lastly, we 
examine the physical causes of the steep slope of the luminos- 
ity scaling relations, concluding that variations of gas content 
with total mass are most likely the dominant reason why these 
are steeper than expected. Our Appendix details the REXCESS 
survey volume calculations and a first attempt at correcting for 
Malmquist bias in the luminosity-mass relation. 

We adopt a ACDM cosmology with Hq = 70 km s"^ Mpc"\ 
Qm = 0.3 and Qa = 0.7, and all uncertainties are quoted at the 
68 per cent confidence level. All logarithmic quantities are given 
to base e, and the quantity L refers to the bolometric [0.01-100 
keV] X-ray luminosity. 

2. Data analysis 

Full details of the sample, including XMM - Newto n observation 
details, can be found in Bohringer et alj (l2007 l). Two of the 
REXCESS clusters, RXC J0956.4-1004 (the Abell 901/902 su- 
percluster) and J2 157.4-0747 (a bimodal cluster), display com- 
plex morphology and are excluded from the present analysis. 



The basic characteristics of the clusters discussed in this paper 
are given in Table [T] 



2.1. Scaling 

In order to estimate cluster quantities consistently, we define 
them in terms of the radius at which the mean mass density 
is 500 times the critical density at the cluster redshifQ. While 
R500 can be estimated from the total mass profile derived un- 
der the assumption of hydrostatic equilibrium (HE), the present 
sample contains clusters in a wide variety of dynamical states 
and consequently the HE a ssumption may n ot be valid in all 
cases (see the discussion in iPratt et"ani2007l) . Instead we esti- 
mate R500 using Yx as a mass proxy. This quantity, defined as 
the product of Mg^soo, the gas mass within R500, and the spec- 
troscopic temperature in the [0.15 - 1] R500 region, is the X-ray 
analogue of the integrated SZ signal Ysz, and has been shown 
to be a low scatter m ass proxy in the numerical simulations of 
iKravtsov et al.l (l2006l) even in the presence of significant dynam- 
ical activity. Recent observational investigations using a variety 
of cluster samples h ave demonstrated that Yx is indeed a low- 
scatter mass proxy (iMaughanI l2007l: lArnaud et^aP l2007h . and 
the theoretical re sults have been ve r ified in independe nt numer- 
ical simulations (iPoole et al.ll2007l: I Yang et al.ll2008h . We esti- 
mate R500 iteratively from the Msoo-Yx relation derived from 
XMM-Newton observations of a sample of 10 nearby morpho- 
logically relaxed local clusters by lArnaud et al.l (l2007l) . viz.. 
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hvoMo, 



(1) 



which was derived using substantially similar methods to those 
described in this paper. The REXCESS gas density profiles from 
which Mgas is derived are discussed in Croston et al. ( 2008). 

Note that there is an ~ 8 per cent normalisation offset of 
the observed relati on when compared to the relation derived by 
iNagai et aP (l2007h from numerical simulations. However, an it- 
erative measurement of R500 from the simulated Msqq-Yx rela- 
tion changes the values of the temperature and luminosity by less 
than 1.5 per cent on average, due to the steep drop of emission 
with radius. Simulations also suggest a ±8 per cent scatter about 
the MsoQ-Yx. Using randomisation assuming a 1.5 per cent rel- 
ative change in the measured quantities due to this scatter, we 
have verified that the slopes and normalisations of the scaling 
laws do not change, and that the maximum change in the scatter 
about the relations is only 7 per cent. 

2.2. Luminosities and temperatures 

Bolometric X-ray luminosities, referred to as L throughout the 
remainder of this paper, were derived for two apertures: (i) the 
entire cluster emission interior to R500 (hereafter Li) and (ii) in 
the [0.15-1] 7?5oo aperture (hereafter L2). We estimated the count 
rates from surface brightness profiles in the [0.3-2] keV band, 
and used the best fitting spectral model estimated in the same 
aperture to convert the count rate to bolometric ([0.01-100] keV) 
luminosity. At 3(T significance, the surface brightness profiles 
are detected out to at least 0.87?5oo for all clusters. For 1 1 clusters, 
we extrapolated the surface brightness profile using a power law 
with a slope measured from that of the data at large radius. As 
can be seen in Table[T] in most cases the need for extrapolation is 



1 M500 = 500p,(z)(4;r/3)/?3^o, where p,(z) 



hHz)3Hl/S7TG and 
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Table 1. Cluster properties. The luminosity is the bolometric [0.01-100] keV luminosity. All quantities are calculated assuming = 0.3, 
Qa = 0.7, and ho = 0.7. 



Cluster 


z 


Ti 


Li 


T2 


L2 


T3 


Yx 


R500 


^det 


cc 


Disturbed 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



RXCJ0003, 
RXCJ0006, 
RXCJ0020, 
RXCJ0049, 
RXCJ0145 
RXCJ0211 
RXCJ0225 
RXCJ0345 
RXCJ0547 
RXCJ0605 
RXCJ0616 
RXCJ0645 
RXCJ0821 
RXCJ0958 
RXCJ1044, 
RXCJ1141 
RXCJ1236 
RXCJ1302, 
RXCJ1311 
RXCJ1516 
RXCJ1516 
RXCJ2014, 
RXCJ2023, 
RXCJ2048, 
RXCJ2129 
RXCJ2149 
RXCJ2157 
RXCJ2217 
RXCJ2218 
RXCJ2234 
RXCJ2319 



8+0203 
0-3443 
.7-2542 
.4-2931 

0- 5300 
4-4017 

1- 2928 
.7-4112 
.6-3152 
.8-3518 
.8-4748 
.4-5413 
.8+0112 
.3-1103 
.5-0704 
.4-1216 
.7-3354 
.8-0230 
.4-0120 
.3+0005 
.5-0056 
.8-2430 

0- 2056 

1- 1750 
8-5048 
.1-3041 
.4-0747 
.7-3543 
.6-3853 
.5-3744 
.6-7313 



0.0924 
0.1147 
0.1410 
0.1084 
0.1168 
0.1008 
0.0604 
0.0603 
0.1483 
0.1392 
0.1164 
0.1644 
0.0822 
0.1669 
0.1342 
0.1195 
0.0796 
0.0847 
0.1832 
0.1181 
0.1198 
0.1538 
0.0564 
0.1475 
0.0796 
0.1184 
0.0579 
0.1486 
0.1411 
0.1510 
0.0984 




876.7 
1059.3 
1045.3 
807.8 
1089.3 
685.0 
693.9 
688.4 
1133.7 
1045.9 
939.2 
1280.0 
755.9 
1077.4 
931.9 
885.2 
753.5 
842.1 
1319.2 
989.9 
927.0 
1155.3 
739.5 
1078.0 
900.6 
886.6 
751.5 
1022.6 
1130.1 
1283.2 
788.7 



0.84 
0.93 
1.07 
0.93 
1.23 
1.33 
0.91 
0.89 
1.32 
1.17 
1.12 
1.28 
0.93 
0.78 
1.09 
1.25 
0.99 
1.22 
1.31 
1.29 
1.37 
1.09 
0.86 
1.48 
0.93 
1.26 
0.97 
1.33 
1.04 
1.15 
1.11 



/ 
/ 
/ 

/ 
/ 



/ 

/ 
/ 

/ 



/ 

/ 

/ 
/ 
/ 

/ 

/ 

/ 



Columns: (1) Cluster name; (2) z: cluster redshift; (3) Ti : spectroscopic temperature of the R < R500 region in keV; (4) Li : luminosity in the R < R500 
region in units of 10"^^ erg s~^ ; (5) T2: spectroscopic temperature in the [0.15-1] R500 region in keV; (6) L2 : luminosity in the [0.15-1] R500 region 
in units of 10"^^ erg s~^; (7) T^: spectroscopic temperature in the [0.15 - 0.75] Rsoo region in keV; (8) Yx in units of 10^^ M© keV; (9) R500 in kpc; 
(10) ratio of the detection radius of the surfa ce b rightness profile at 3cr significance, T^det, to R500; (11) systems classified as cool cores on the basi s 
of central density vs. cooling time (see Sect. l2.3l ); (12) systems classified as disturbed on the basis of the centre shift parameter (w) (see Sect. l2.3l ). 



minimal. IVEeasured luminosities for both apertures are given in 
Table [T] Errors take into account statistical factors, uncertainties 
in R500 and extrapolation uncertainties. These were estimated 
from IMonte Carlo realisations in which the above procedure, 
including extrapolation, was repeated for 100 surface brightness 
profiles, the profiles and R500 values each being randomised ac- 
cording to the observed uncertainties. A PSF correction was im- 
plemented by using the gas density profile derived from reg- 
ularised dep rojection of the surface brightness as described in 
ICroston eTa l. (2006); the correction was obtained from the ratio 
of the observed to PSF-corrected count rates in each aperture. 
The correction is negligible in the full aperture but can be up to 
13 per cent for strong cooling core systems in the core-excluded 
aperture. 



Spectroscopic temperatures were measured in the [0.15 - 
0.75] /?5oo apertur43 (referred to hereaf ter as 7^3) from iter - 
ation about the M500 - Yx relation of lArnaud et al.l (l2007l) . 
Temperatures in the full aperture and in the [0.15 - 1] R500 aper- 
ture (hereafter Ti and T2, respectively) were then calculated by 
re-extraction of spectra in the relevant regions. In all cases the 
spectra were fitted with a IMeKaL model with an absorption fixed 
at the HI value (excepting RXC J2014.8-2430, which was found 
to have a significantly higher absorption than that indicated from 
the HI value). The three EPIC cameras were fitted simultane- 
ously in the [0.3-10] keV band, with the regions around the in- 
strumental lines (1.4-1.6 keV for all cameras and 7.45-9 keV for 
the pn camera) excluded from the fit. Temperatures for all three 
apertures are listed in Table [T] 



We have also calculated luminosities in the [0.1 - 2.4] keV 
and [0.5-2] keV bands for ease of comparison with previous 
soft X-ray survey results from ROSAT. Luminosities and their 
associated scaling relations are given in Appendix|Bl 



^ Note that this aperture is more appropriate for comparison with dis- 
tant clusters, which often have poor signal to noise in the outer regions. 
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Table 2. Observed bolometric X-ray luminosity scaling relations. For each set of observables (L,A), we fitted a power law relation of the form 
hizTL = C (A/AoY, with Aq = 5 keV, 2 x 10^^ keV and 2 x 10^^ M©, and n = -I, -9/5 and -7/3 for T, Yx and M, respectively. Results are 
given for the BCES (Y|X) and BCES orthogonal fitting methods (see Section [241 ). The intrinsic natural logarithmic scatter about the best fitting 
relation in the In-ln plane is also given for each fit. ^ Since M is derived from Yx, the values of the scatter in the L - M relation are identical to 
those for the L- Yx relation. ^ Corrected for Malmquist bias (see Appendix 151). 



Relation Fitting Method 







BCES (Y|X) 




BCES Orthogonal 




C(10^^ergs-0 


a 


<^ln L,intrinsic 


C(10^^ ergs-^ 


(X 


L,intrinsic 


All 






R < R500 






L1-T3 
Li-Yx 
Li-My 
Li-My MB^ 


6.07 + 0.58 
5.62 + 0.46 
5.20 + 0.36 
1.81+0.13 
1.45 + 0.12 


2.70 + 0.24 
2.88 + 0.23 
0.99 + 0.05 
1.81+0.10 
1.90 + 0.11 


U.DDj ± U. i ID 

0.525 + 0.097 
384 + 060 

\J»^KJ— — !— V7.V7V7V7 

a 
a 


7 1 a _i_ 1 H'^ 

6.27 + 0.67 
5 35 + 38 
1.74 + 0.13 
1.38 + 0.12 


Q + '^9 

D.DJ X \J.DZ^ 

3A2 ± 0.27 
1.04 + 0.06 
1.96 + 0.11 
9 OS + n 1 ^ 

Z.Oo It 0. i J 


yj. / DD VJ.UJ 

0.560 + 0.115 
0.383 + 0.061 

a 
a 


Cool core 














Li-Ti 
Li-Yx 


11.15 + 2.42 
7.71+0.58 


2.71 + 0.48 
1.04 + 0.07 


0.432 + 0.108 
0.234 + 0.103 


12.79 + 3.80 
7.84 + 0.65 


3.15 + 0.63 
1.06 + 0.09 


0.479 + 0.135 
0.236 + 0.107 


Non-cool core 














Li-Ti 
Li-Yx 


4.78 + 0.29 
4.27 + 0.20 


2.89 + 0.21 
0.96 + 0.05 


0.267 + 0.058 
0.214 + 0.035 


4.97 + 0.29 
4.32 + 0.20 


3.06 + 0.19 
0.98 + 0.05 


0.285 + 0.068 
0.214 + 0.036 


Disturbed 














Li-Yx 


4.18 + 0.59 
3.72 + 0.27 


2.49 + 0.56 
0.92 + 0.09 


0.497 + 0.215 
0.245 + 0.120 


5.43 + 2.74 
3.85 + 0.32 


3.19 + 0.78 
0.96 + 0.08 


0.646 + 0.346 
0.249 + 0.123 


Regular 














Li-Yx 


7.26 + 0.86 
6.15 + 0.42 


2.62 + 0.21 
0.97 + 0.05 


0.578 + 0.118 
0.302 + 0.058 


7.97 + 1.28 
6.21 + 0.44 


3.13 + 0.33 
1.00 + 0.05 


0.634 + 0.142 
0.303 + 0.059 



All 








1^2-7^2 


3.89 + 0.18 


2.78 + 0.13 


0.269 


^2-7^3 


3.31+0.16 


2.84 + 0.17 


0.331 


L2-YX 


3.05 + 0.07 


0.97 + 0.03 


0.156 


L2-MY 


1.09 + 0.05 


1.77 + 0.05 


a 


Cool core 








L2-T2 


4.31+0.42 


2.58 + 0.23 


0.242 


L2-YX 


3.36 + 0.16 


0.96 + 0.04 


0.144 


Non-cool core 








L2-T2 


3.74 + 0.21 


2.89 + 0.18 


0.231 


L2-Yx 


2.91 + 0.06 


0.97 + 0.03 


0.114 


Disturbed 








L2-T2 


3.58 + 0.41 


2.88 + 0.37 


0.295 


L2-YX 


2.77 ± 0.07 


0.99 ± 0.04 


0.111 


Regular 








^2-7^2 


4.13 + 0.21 


2.68 + 0.11 


0.225 


L2-Yx 


3.24 + 0.08 


0.94 + 0.02 


0.115 



0.15 <R<R5oo 



+ 0.055 
+ 0.068 
+ 0.038 


4.06 + 0.22 
3.48 + 0.21 
3.06 + 0.07 
1.08 + 0.04 


2.94 + 0.15 
3.07 + 0.18 
0.98 + 0.03 
1.80 + 0.05 


0.279 + 0.059 
0.346 + 0.075 
0.156 + 0.038 

a 


+ 0.110 
+ 0.098 


4.46 + 0.56 
3.38 + 0.17 


2.70 + 0.26 
0.97 + 0.05 


0.247 + 0.113 
0.145 + 0.098 


+ 0.035 
+ 0.027 


3.88 + 0.22 
2.92 + 0.06 


3.02 + 0.19 
0.98 + 0.03 


0.237 + 0.039 
0.114 + 0.027 


+ 0.080 
± 0.096 


4.00 + 0.73 
2.79 ± 0.08 


3.18 + 0.38 
0.99 ± 0.04 


0.312 + 0.090 
0.111 + 0.096 


+ 0.070 
+ 0.045 


4.20 + 0.23 
3.24 + 0.08 


2.76 + 0.11 
0.94 + 0.02 


0.231 +0.075 
0.115 + 0.045 



Li/Ti : luminosity/temperature interior to/?5oo; L2/T2: lumi nosity/temperature in the [0.15-1] R500 aperture; T3: temperature in the [0.15-0.75] R500 
aperture; My: mass measured from the M^qq-Yx relation of Arnaud et al.l (l2007h . 



2.3. Subsamples 



2.3.1. Cooling time classification 



We further subdivide the sample to elucidate the eflTects of cool 
cores and merger-related phenomena on the scaling relations. 



In an approach similar to that used by lO'Hara et 10(12006'), we 
use gas d [ensity and cooling ti me profiles to classify cooling core 
systems. ICroston eFaP (l2008h describe the gas density and cool- 
ing time profiles of the present sample, which are fully depro- 
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0.001 0.010 0.100 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.01 0.10 1.00 

h(z)"^ n,_o (cm"^) log^o (w> ([0.15-1] R^oo) Radius (R500) 



Fig. 1. Definition of cluster subsamples. Left: Central cooling time vs. central gas density, Uefi. The dotted line delineates the threshold we use to 
define cool core systems: h(z)~^nefi > 4 x 10~^ cm~^, ^cooi < 10^ years. Centre: Histogram of centre shift parameter (w>, evaluated in the [0.15- 
l]i?50o aperture. Clusters with {w) > 0.01 /?5oo are classified as morphologically disturbed. Right: Emission measure profiles of the R EXCESS 
sample, scaled according to the standard dependence on temperature and expected evolution with redshift. Systems classified as cool core and as 
morphologically disturbed are indicated (see Sect. 12.31 1. 



jected an d PSF-corrected usin g the non-parametric method de- 
scribed in Crosto n et al ] (l2006l) . We estimate the central gas den- 
sity Hefi from a p model fit to the deconvolved, deprojected gas 
density profiles interior to 0.037?5oo. Figure [T] shows the central 
cooling time versus ^^^o. which exhibits a strong correlation, as 
expected since the cooling time is derived from the gas density. 
We classify clusters according to their central gas density ^^,0, 
such that those with h{z)~^ rie^o > 4 x 10"^ cm"^ (equivalent to 
those with a central EM value EMc^n > 20 x 10"^ cm"^ Mpc in 
the right hand panel of Fig. [T]) are defined as cool core systems; 
10/31 clusters are classified as such. Figure [T] shows that these 
systems have central cooling times ^cooi,o < 10^ years. 

2.3.2. Morphological classification 

The samp le also contains clusters in a wide vari ety of dynami- 
cal states (iBohringer et al .120071: [Pratt et al.l2'007h . To investigate 
the eflfect of dynamical state on the relations, we have calculated 
values of centroid shift (w), defined as the standard deviation of 
the projected separations between the X-ray peak and the cen- 
troid at each radius in the [0.1-1] R500 region: 



<w) = 



i^^(A.-(A))^ 



1/2 



1 

^500' 



(2) 



where A/ is the projected distance between the X-ray peak and 
centroid in the / th aperture. 

Introduced by 'Mohr et aP (Il993h . this quantity was found 
to be the most sensitive i ndicator of dynam ical activity in the 
numerical simulations of Poole et aD (l2007h . We calculate cen- 
troids in circular apertures of radii ^ x 0.1 x R500 with n = 
2, 3 ... 10, excluding the central regions to avoid biases associ- 
ated with enhanced emission from cool cores (although exclu- 
sion of the central region does not have a significant eflfect on 
the results). The centroid shift (w) is then defined as the stan- 
dard deviation of the projected separations between the X-ray 
surface brightness peak and the centroid in units of Rsoo- A forth- 
coming paper will discuss these results in more detail. For the 
current analysis, the distribution of (w) for the present sample 
is shown in the central panel of Figure [T] We classify clusters 
with (w) > 0.01 7?5oo as morphologically disturbed, and clusters 



with (w) < 0.01 R500 as morphologically regular. In total, 12/31 
clusters are defined as morphologically disturbed. 

The diflTerent subsample classifications are indicated in 
Table [T] and illustrated in Figure [T] In general, the presence of 
a cool core is anti-correlated with indications for morphological 
disturbance. However, two clusters possess both a cool core and 
display evidence for morphological disturbance (RXCJ1302.8 
-0230 and RXCJ2319.6 -7313). A gallery of the cool core 
and non-cool core systems, sorted by (w), can be found in 
Figures IaTT] and IA.21 respectively, in Appendix lAl 

2.4. Fitting procedure 

For each set of observables (B,A), we fitted a power law rela- 
tion of the form h(zTB = C(A/Ao)^, where h(z) is the Hubble 
constant normalised to its present day value and n was fixed to 
the expected scaling with z. The fit was undertaken using linear 
regression in the log-log plane, taking the uncertainties in both 
variables into account. Assuming a linear relation of the form 
Y = aX-\-b, and a sample of N data points (7/, X/) with errors cry. 
and (TXi , the raw scatter was estimated using the error weighted 
orthogonal distances to the regression line: 



1 ^ 

i-i 



where 



Wi 



l/cr 



■ aXi - by 



and 



2 22 
o-y, + a (T^.. 



(3) 



(4) 



The intrinsic scatter was computed from the quadratic diflTerence 
between the raw scatter and that expected from the statistical 
uncertainties. 

As Figures [2]- [5] show, the uncertainties in the present data set 
are entirely negligible compared to the intrinsic scatter, so that 
error weighting of individual data points will have no effect on 
the resulti ng fits. In the following w e use the BCES regression 
method ( Akritas & Bershadylll996b . which takes into account 
measurement errors in both coordinates and intrinsic scatter in 
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the data and is widely used in astronomical regression, giving 
results that may easily be compared with other data sets fitted 
using the same method. 

It is well-known that diff'erent regression methods give dif- 
ferent slopes even at the po pulation level (e.g.Jisobe et al.llT99Ql: 
lAkritas & Bershadylll996h . It is therefore of paramount impor- 
tance to choose the regression method best suited to the data in 
hand. With the present data set, there is no easy answer to the 
question of which quantity to treat as the dependent variable and 
which to treat as the independent variable. In cosmological and 
theoretical applications, the mass of a cluster is its most funda- 
mental property. Given the tight mass-temperature relation (e.g., 
lArn aud et al. 2005), it is reasonable to assume that T is closely 
coupled to the mass. However, as will be seen below, there is a 
large intrinsic scatter in L, presumably due to baryon physics. 
One possible minimisation method would thus treat L as the de- 
pendent variable. A second possible minimisation method would 
be to assume that both variables are quasi-independent, and to 
treat them symmetrically. 

In the following, we thus give the results from the BCES 
(Y|X) fitting method, which minimises the residuals in L, and 
from the BCES orthogonal fitting method, which minimises the 
squared orthogonal distances. In the case of maximum scatter 
(the raw, uncorrected Li-Ti relation), the BCES (Y|X) method 
typically gives slightly shallower slopes than the orthogonal 
BCES metho(£|. As the scatter decreases, the various regression 
methods give results which agree very well within their Icr un- 
certainties (Table [21). Uncertainties on all fit parameters and as- 
sociated scatter are determined from 10 000 bootstrap resamples 
of the data. Since measurement errors are at the 1 - 3 per cent 
level, we give only estimates of the intrinsic dispersion about the 
best fitting relations. 

3. Results 

3.1. Scaled emission measure profiles 

The emission measure (EM) was calculated from the surface 
brightness profiles extracted in the [0.3-2] keV band via: 

j7A4r ^ 47r(l+z)^^(^(x)) 

EM(r) = — ; r = dA(z)0, (5) 

^(T, z) 

where S(6) is the surface brightness, Ja(^) is the angular dis- 
tance at redshift z, and e(T, z) is the emissivity, which has been 
calculated t aking into account ab sorption and the instrument re- 
sponse (e.g. [Neumann & Arnaudl ll999). We then scaled the EM 
profiles according to their expected evolution with redshift and 
dependence on temperature, EM oc h(z)~^T~^^^, shown in the 
right hand panel of Figured] The behaviour of the scaled profiles 
is ver y similar to that seen in the gas density profiles discussed 
in ICroston et al.l (l2008b : outside the central regions, the disper- 
sion rapidly decreases and the profiles begin to show indications 
of similarity. The relative dispersion in scaled profiles shows a 
broad minimum of cr/{EM(r)} ~ 0.35 from 0.2 - 0.97?5oo, with 
a maximum of 1.56 in the central regions and a minimum of 
0. 32 at 0.5/?5oo. The lat t er is s omewhat smaller than that found 
by [Neumann & ArnaudI (Il999h . who used a relation taken from 
numerical simulations to calculate R500. 

The two subsamples form distinct classes in the plot. The 
cool core systems, unsurprisingly, show very strong central 

^ The BCES (Y|X) method gives precisely the same results as 
the modifi ed weighted least squares (WLSS) method described in 
lPrattetaLl ( [2006) . 



emission and also appear to scale somewhat more tightly in the 
outer regions. In contrast, all of the clusters with the lowest cen- 
tral emission measure are classed as disturbed. 

3.2. TheL-T relation 

In the left-hand panel of Figure [2] we show the Li-Ti rela- 
tion obtained with Ti and bolometric Li derived from all emis- 
sion interior to R500 (i.e., equivalent to a raw, uncorrected re- 
lation). In many cases the errors are smaller than the points, 
a testament to the exceptional quality of the data. The best 
fitting power law relations derived from the BCES fits are 
overplotted; fits are Hsted in Table [2l The BCES (Y|X) slope, 
2.70 ± 0.23, is consistent with previou s determinations such as 
those of'Markevitch"(1998', 2.64+0.16), Arnaud & Evrard (19991 
2 .88 + 0.15), Allen & Fabian (1998, 2.9 + 0.3) and Novicki et al] 
(120021 2.82 + 0.32). The slope derived from the BCES orthogo- 
nal fit, 3.35 + 0.32 is somewhat steeper, although only at slightly 
more than 2cr, a consequence of th e very large scatter in the data . 
The relations of Arnaud & Evrard' (T999) and Markevitch[ ([1998[) 
are also plotted in the Figure: their normalisations are notably 
lower that that found in the present work, due to their being a 
non-cool core cluster sample and cool core-corrected sample re- 
spectively. 

The central panel of Fig. [2] shows the histogram of the log 
space residuals from the best fitting relations. The logarithmic 
scatter about the Li-Ti relation is cr\^L ~ 0.7 + 0.1 in both cases, 
and is dominated by the intrinsic component, the statistical scat- 
ter being negligible. Although the residuals present a clear skew 
towards higher luminosity systems, the KS probability that the 
residuals are drawn from a Gaussian distribution are 0.15 and 
0.10 for the BCES (Y|X) and orthogonal fits, respectively. This 
result does not strongly exclude the Gaussian hypothesis, un- 
derlining the need for a la rger sample to bett er understand the 
scatter (although note that [Novicki et al.[ '2002 find that the log 
space residuals of a larger sample are consistent with a Gaussian 
distribution with a similar cr to that found for the present sam- 
pie). 

It is interesting to investigate the factors driving the very 
large scatter about the Li-Ti relation. The right hand panel of 
Figure [2] shows the log space deviations from the best fitting 
Li-Ti relations for the cooling core sample (blue stars), and 
for the morphologically disturbed subsample (red squares). In 
both cases, the subsamples clearly populate diff'erent regions of 
the residual space: cool core systems are preferentially located 
above the main relation, while morphologically disturbed sys- 
tems lie below it. The best fitting power law relations to the indi- 
vidual subsamples are listed in Table[2l Cool core clusters have a 
statistically identical slope to that of the non-cool core systems, 
and to that of the sample as a whole. The cool core subsample 
has a higher normalisation than the non-cool core subsample, 
significant at the > 2cr level, suggesting that the primary eff'ect 
of a cooling core is to move a given system orthogonally from 
the standard relation. However, the logarithmic scatter about the 
cooling-core only relation (cr\^L = 0-48 + 0.13 for the BCES 
orthogonal fit) is higher than that about the non-cool core re- 
lation at slightly more than Icr, reflecting the wide variety of 
cooling core strengths in the present sample. The logarithmic 
scatter about the non-cool core relation ((XinL = 0.29 + 0.07, 
or ^logio ~ 0.15), is in good agreement with that found by 
'Arnaud & Evrard ([1999|) . from a sample which contained only 
non-cool core systems. 

A similar trend is seen when the clusters are divided accord- 
ing to the morphology parameter. Firstly, it is clear from Figure[2| 
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Fig. 2. Left: Li-Ti relation for the R EXCESS sample (quantities derived from all emission interior to Rsqq). The error bars are smaller than the 
points in many cases. The best fitting power law relation derived from the BCES (Y|X) (red line) and BCES orthogonal (blue line) are overplotted; 
the shaded r egion corresponds to the Icr uncertainty on the latter. The dashed and dot-dashed lines are the relations of Arnaud & Evrard ( 1999) and 
iMarkevitchI ([l998), respectively. Centre: Histogram of the log space residuals from the best fitting L-T relation, derived from each fitting method 
as indicated. Right: Log space residuals for both fitting methods as indicated. Cooling core clusters (blue stars) and morphologically disturbed 
clusters (red squares) occupy two distinct regions in the plot in both cases. 
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Fig. 3. Left: L2-T2 relation for the REXCESS sample (quantities derived from emission in the 0.15/?5oo < R < R500 aperture). The best fitting 
power law relation derived from the BCES orthogonal fitting method is overplotted as a solid line (the BCES ( Y|X) results are ve r y sim ilar); 
the shaded region corresponds to the Icr uncertainty on the fit. The dashed and dot-dashed lines are the relations of lArnaud & EvrardI 41999) and 
iMar kevitch (1998), respectively. Centre: Histogram of the log space residuals from the best fitting L-T relation, derived from the BCES orthogonal 
fit method. The solid curve is a Gaussian with crx^i - 0.28, corresponding to the scatter about the best fitting relation. Right: Log space residuals. 
Cooling core clusters (blue stars) and morphologically disturbed clusters (red squares) are less obviously segregated once the central region is 
excised. 



that the disturbed clusters preferentially populate the lower enve- 
lope of the Li-Ti relation. The slopes of the relations are statis- 
tically identical for both subsamples, and in agreement with that 
of the entire sample, but the normalisation of the relaxed sam- 
ple is higher at the Icr level than that of the unrelaxed sample. 
This is partly due to the predominance of cool core systems in 
the relaxed subsample, although disturbed cool core systems do 
exist. The logarithmic scatter about the relations is very similar, 
at CTinL ~ 0.65, although they are not well constrained. 

The clear segregation of the cooling core clusters from the 
rest of the population (Fig. [2]), together with the small relative 
segregation of dynamically disturbed systems and the structural 
similarity at large radius (Fig. [T] right hand panel), suggest that 
simply excluding the central region should tighten the luminos- 
ity scaling relations. Figure [3] shows the L2-T2 relation derived 



from emission excluding the core region, where both the lumi- 
nosity and temperature are estimated in the 0.15 R500 < R < R500 
aperture; best fitting slopes and normalisations are given in 
TableEl 

The BCES (Y|X) slope, 2.78 ± 0.13 is similar to the relation 
for that derived from all emission interior to R500, and the slope 
of the BCES orthogonal fit, 2.94 ± 0.15 is slightly steeper but 
in good agreement within the Icr uncertainties, as is the normal- 
isation. This relation is in excellent agreement, both in terms 
of slope and normalisation with those of IMarkevitch, (1998<) 
an d iArnaud & EvrardI (1 19991) . which are overplotted in the same 
Figure. 

However, the logarithmic intrinsic scatter, crinL = 0.27 + 0.06 
is smaller, as expected, by a factor of two. The central panel of 
Figure [3] shows the histogram of the log space residuals from 
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Fig. 4. Left. Li-Yx relation for the REXCESS sample, with luminosity derived from all emission interior to Rsqq. The error bars are smaller than 
the points in many cases. The best fitting power law relation derived from the BCES orthogonal fitting method is overplotted as a solid line (the 
BCES (Y|X) results are very similar); the shaded region corresponds to the Icr uncertainty on the fit. The dashed line is the fit derived by MaughaiJ 
(l2007h from observations of 115 galaxy clusters in the Chandra archive. The agreement is excellent. Centre: Histogram of the log space residuals 
about the best fitting orthogonal BCES relation. Right: Log space residuals of the diff'erent subsets. Blue stars: cooling core clusters; red squares: 
morphologically disturbed systems. 
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Fig. 5. Left: L2-Yx relation for the REXCESS sample, with luminosity derived from emission in the 0.15 /?5oo < R < R500 aperture. The best 
fitting power law relation derived from the BCES orthogonal fitting method, which takes into account errors in both coordinates and intrinsic 
scatter in the data, is overplotted as a solid line (the BCES (Y|X) results are very similar); the shaded region corresponds to the Icr uncertainty on 
the fit. The dashed line is the fit derived by Maughan ( 2007) from observations of 115 galaxy clusters in the Chandra archive. The agreement is 
again excellent. Centre: Histogram of the log space residuals about the best fitting orthogonal BCES relation. The solid curve is a Gaussian with 
cfinL - 0.16, corresponding to the scatter about the xoXdiiion.Right: Log space residuals of the diff'erent subsets. Blue stars: cooling core clusters, 
red squares: morphologically disturbed systems. 



the best fitting BCES orthogonal relation. The overplotted curve 
is a Gaussian with raw cr\^L = 0.27, corresponding to the scat- 
ter about the relation. The KS test probability the residuals are 
drawn from a Gaussian distribution is 0.54 for the BCES orthog- 
onal fit. 

The scatter is clearly reduced on exclusion of the core re- 
gions. The relative eflfect of the change in luminosity and tem- 
perature in the reduction of scatter can be estimated simply by 
comparing the values estimated in the two apertures. We find 
{T1/T2) = 1.02 ± 0.07 and {L1/L2) = 1.62 + 0.31 for the full 
sample, indicating that the change in temperature is a very mi- 
nor effect compared to the change in luminosity. Unsurprisingly 
however, the change in temperature is negative for cool core sys- 
tems ({T1/T2} = 0.96 ± 0.07), while it is positive for non-cool 
core objects ({T1/T2} = 1.05 ± 0.05). 



Table [2] also lists the fits to the diflTerent subsamples. Scatter 
decreases markedly (by approximately a factor of two) for the 
cool core subsample, as expected, but it also decreases some- 
what (-15 per cent) for the non cool core subsample. Cool core 
clusters still tend to be found towards the upper envelope of the 
distribution, which is reflected in their slightly higher normalisa- 
tion compared to non-cool core systems (although this is not sig- 
nificant). The slopes are stable however, and in agreement with 
those found for the relation derived from all emission interior 
to 7?5oo. In common with the full emission sample, morpholog- 
ically disturbed clusters tend to describe the lower envelope of 
the distribution, having a slightly lower normalisation than the 
full sample, although this is not significant, and a similar slope. 



G.W. Pratt et al.: X-ray luminosity scaling relations of the R EXCESS 



9 




M^oo (h^"o M J M^oo (h^"o M^oo (h^"o Mj 



Fig. 6. L - M500 relation for the R EXCESS sample, with the mass estimated from the Yx - M relation of lArnaud et alJ (l2007l) . Left: Relation 
for all emission interior to R500. Centre: Relation for emission in the [0.15 - 1] R^qq aperture. The best fit t ing po wer law relation derived from the 
orthogonal BCES fit method is overplotted as a solid line. The dashed line is the fit derived by'Maughan' (l2007l) from observations of 115 galaxy 
clusters in the Chandra archive. Right: Relation corrected for Malmquist bias as described in Appendix IbI 



3.3. The L-Yx relation 

Yx is an interesting quantity because simulations suggest that 
deviations in Mgas,5oo and T for a given system are anti- 
correlated with respect to the self- similar expectations, leading 
to a reduction in scatter (although thus far this is empirically 
untested). That Yx is the X-ray analogue of the integrated SZ 
Comptonisation parameter Ysz makes the calibration of its rela- 
tionship with the X-ray luminosity of prime importance for the 
interpretation of data from the upcoming all-sky surveys from 
the Planck and eROSITA satellites. 

The Li-Yx relation, where the luminosity is derived from 
all emission interior to R500, is shown in Figure (H the best fit- 
ting power law values are given in Table [2l Because of the 
smaller scatter in these data, both BCES fitting methods give 
consistent results. Our relation is in good agreement with that of 
the Ch andra archive study of 115 galaxy clusters by 'Maughan' 
(l2007h : the slope B = 1.10 ± 0.04 is consistent with our BCES 
orthogonal value, a = 1.04 ± 0.06, and the normalisation at 
Yx = 2 X 10^"^ MokeV is only 14 per cent lower than ours. The 
intrinsic logarithmic scatter is (Tin l = 0.38 + 0.06, considerably 
less than about the L-T relation (note that since Yx is calculated 
using the temperature estimated in the [0.15 - 1] R500 aperture, 
this will tend to damp scatter somewhat). However, the residual 
histogram about the best fitting orthogonal BCES relation, plot- 
ted in the central panel of the same Figure, has a KS probability 
of only 0.09 of being compatible with a Gaussian distribution. 

The residual distribution for the diflferent subsamples mirrors 
that of the Li-Ti relation: cool core clusters lie preferentially 
above and morphologically disturbed systems lie preferentially 
below. This fact is reflected in the diflferent normalisations found 
when fitting the diflTerent subsamples: the cool core systems have 
the highest normalisation and the disturbed systems the lowest. 
However, the slope of the relation, when fitted to diflTerent sub- 
samples, is remarkably stable at 0.96-1.06, and the slopes of all 
subsamples are statistically indistinguishable. 

Figure [5] shows the L2-Yx relation, determined with the core 
emission excluded. Once again there is excellent agreement in 
both slope and n ormalisation between our relation and that of 
iMaughan ('2007). The relation is very tight: the intrinsic loga- 
rithmic scatter is only (XinL = 0.16 + 0.04, and the KS probability 
that the distribution of residuals is compatible with a Gaussian is 
0.93. For the diflTerent subsamples, Table[2l shows that the slopes 



are remarkably similar, ranging from 0.94 to 0.99, and the power 
law normalisations for the best fitting models are segregated in a 
similar manner to the luminosity temperature relation, although 
with much reduced significance. 



3.4. The L - M500 relation 

It is interesting to make a first examination of the slope and nor- 
malisation of the L - M500 relation for the present sample. Since 
we do not have independent meas ures o f the mass, we use the 
^500 - Yx relation of Arnau d et aP (|2007l) to estimate the masses 
of the clusters in the sample. For the purposes of this initial in- 
vestigation, we ignore the impact of the intrinsic scatter about 
the M500 - Yx relation because it is at present not sufl&ciently 
well quantified; X-ray calibrations are necessarily available for 
relaxed cluster samples only, and weak lensing calibrations are 
at present lacking sufl&cient dynamic range in mass. The present 
approach allows us to verify the slope and normalisation of the 
relation under the given assumptions, to check the coherence of 
the slopes, and to compare with previous work using similar ap- 
proaches. 

The measured L - M relations are summarised in Table [2] 
and the relations obtained for bolometric L measured in both 
apertures are plotted in Figure [6] The slopes of the relations, 
~ 1 .8, are consistent with the L-T and M500 - Yx relations, as ex- 
pected. Comparing our measurements of the slope and normal- 
isation with those of Maughan' (^2007), we find excellent agree- 
ment in slope for the relation derived from all emission inte- 
rior to /?5oo, athough our normalisation is somewhat higher (by 
< 20 per cen t). When the core emission is excluded, the slope 
of lMaughanl 's relation (1.63 + 0.08) is somewhat shallower than 
our BCES orthogonal measurement (1.80 + 0.05), at the ~ 2cr 
level. However, the normalisations are in excellent agreement. 

In the right hand panel of Figure[6]we compare the raw L-M 
relation with that corrected for the eflTect of Malmquist bias. The 
correction procedure, and the relations for the [0.1 - 2.4] keV 
and [0.5 - 2] keV ba nds, plus comparison with the results of 
IVikhlinin et aT (2008), are given in Appendix |B] The correction 
has the eflTect of steepening the relation slightly due to the under- 
representation of low-luminosity clusters on the REXCESS sam- 
ple. 
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Table 3. Best fitting parameters for the three parameter scaling re- 
lation fits. Data were fitted with a power law of the form h{zfL - 
C {AlAoYinef, with Aq = 5 keV, 2 x 10^^ keV and 2 x 10^^ M©, 
and n--\, -9/5 and -7/3 for T, Yx and M, respectively. 



Relation C a P crinL,intrinsic 

Li-Ti-Hefi 27.45 + 1.45 2.61 + 0.36 0.36 ± 0.10 0.47 ± 0.04 
Li-Yx-Uefi 13.90 + 1.13 0.99 + 0.04 0.26 + 0.03 0.22 ± 0.02 
Li-M-Hefi 4.84 ±1.14 1.82 ±0.07 0.26 ± 0.03 ... 

Li/Ti: luminosity/temperature interior to Rsqq. 

The scatter is, by definition, identical to that about the L - 
Yx relation, and is i n excellent agreement with that found by 
IVikhlinin et alJ (l2QQ8t) from a similar analysis of a larger flux- 
limited sample of nearby clusters. Note that if L oc M^, then 
a first order estimate of the scatter in mass is (XinM ~ o-ini/y- 
However, this will only be true if (XinL is measured at fixed M for 
a complete sample. Using the measurement of (XinL at fixed T or 
Yx introduces covariance of T and Yx into the relations, which 
would modify the first order scatter estimate. Nevertheless, the 
scale of this first order estimate of the scatter in mass is (XinM ~ 
0.20 - 0.37 for the fufl aperture and only crinM ~ 0.09 - 0.16 for 
core extracted quantities. 

3.5. Relations including a third parameter 

The presence of a cool core is clearly the factor which con- 
tributes most to the scattering of a given cluster about the best fit- 
ting relation. Figure[T] shows that the central density ^ ^.p is a very 
reliabl e indicator of cool core strength. Following Q' Hara et al.l 
(l2006l) . it thus follows that n^^o may be taken into account as a 
third parameter in the scaling relations. 
Fitting a scaling relation of the form: 

KzTL = C(A/Aonnef (6) 

where h(z) is the Hubble constant normalised to its present day 
value and n was fixed to the expected scaling with z, and solving 
for a,l3 and the normalisation C, allows us to investigate the in- 
fluence of central density ne,o on the scaling relations. For each 
relation the fit was undertaken using standard linear regression 
in the log-log plane. We determine the best fitting values and as- 
sociated Icr uncertainties via 1000 bootstrap resamplings of the 
observed data set, and the raw scatter was estimated using the 
error weighted orthogonal distances to the regression line. 

The resulting best fitting relations are summarised in Tabled 
Dependencies on the main scaling parameter (T, Yx and M) 
are similar to those derived for the two-parameter fits to core- 
excluded quantities, as expected. The scatter is comparable to 
that derived from a two-parameter fit to core-excluded quanti- 
ties for all relations. 

Thus the technique of using the central gas density appears to 
be a promising method for reducing scatter about the luminosity 
scaling relations. 

4. Discussion 

4.1. Scatter about the relations, and correcting for it 

The R EXCESS data have allowed us to investigate the sources 
and magnitude of the scatter about the various relations using 
a data set which should be representative of any X-ray selected 
sample of clusters. 



For all emission interior to R500, morphologically disturbed 
systems tend to lie below the best fitting relation to the entire 
whole sample; however, this is mainly due to the eflTect of cool 
core systems in the full sample, which tend to increase its nor- 
malisation. A fairer test is to compare the subsamples when the 
core emission is excluded: in this case, the normalisations dif- 
fer in all cases by less than l.Scr, although morphologically dis- 
turbe d systems sti l l have the lowest normalisation. In common 
with Q'Hara et a n (l2006h . we find that clusters with greater mor- 
phological substructure do not exhibit more scatter about scal- 
ing relations than clusters with less substructure. This is mostly 
a consequence of the fact that morphologically regular systems 
contain a preponderance of cool core clusters. This result would 
also suggest that the main eflTect of merging is to move systems 
along the relation rather than orthogonal to it. In this context, 
we note that numerical simulations predict quasi- simultaneous 
boosting in temperature and luminosity at certain epochs after a 
merging event, which would indeed tend to move clusters along 
the relation. 

In common with most previous investigations, we find that 
the vast majority of the scatter in all relations is due to the pres- 
ence of cool cores, which lie systematically above the best fitting 
relation with an oflTset that appears to be related to the strength of 
the cool core. A fit to the cool core systems diflfers from a fit to 
the whole sample only by a normalisation factor. However, the 
scatter about the best fitting cool core subset relation is nearly 
twice that about the non cool core subset relation, reflecting the 
diflTerent cool core strengths. 

Excluding the central emission leads to a significant reduc- 
tion of the scatter about all relations. For example, the scatter 
about the L-T relation decreases by a factor of two on exclusion 
of the core, and the reduction in scatter is similar for the L - Yx 
relation. Correcting for the presence of a cool core by assuming 
a power law dependence of central density ne,o with luminosity 
aflfords an alternative method to reduce scatter. The reduction in 
scatter obtained by the use of ne,o is of the same order as that ob- 
tained from simple core exclusion. We note that core exclusion 
may be difi&cult in the case of distant clusters or those detected 
with very low signal to noise, and in these circumstances it may 
be preferable to use the central density or surface brightness to 
reduce scatter about the scaling relations. 

4.2. Slope of the relations 

The X-ray luniinosit y of a cluster can be written 
(lArnaud&Evrardlll999l) : 

L(T) = fl,(T)[M(T)A(T)]Q(T) (7) 

where A(T) is the cooling function. Q(T), introduced by 
lArnaud & EvrardI (119991 is equal to <Pgas)/<Pgas)^, with the an- 
gle brackets denoting an intrinsic volume average. Q(T) is thus a 
dimensionless structure factor which depends only on the spatial 
distribution of the gas density (e.g., dumpiness at small scale, 
shape at large scale, etc). With the set of additional assumptions 
(i) pure bremsstrahlung emission [A(r) oc T^^^]; (ii) virial equi- 
librium [M oc r^/^]; identical internal cluster structure [Q(T) = 
Ci]; constant gas mass fraction [fgas(T) = C2], we arrive at the 
standard self-similar expectation for the luminosity-temperature 
relation, LocT^. The self- similar L- relation can be obtained 
from combination of the gas mass-luminosity and luminosity- 
temperature relations to give L oc Y^^^ . Combining the self- 
similar mass-temperature and luminosity-temperature relations 
leads to a dependence of luminosity with mass of L oc M^^^ . 
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T < R^oo (keV) 
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Fig. 7. Structure factor 

2(7^)|bin = <Pgas>/<Pgas>^ es- 
timated from the gas density 
profiles, versus temperature. 
Left: Quantities estimated from 
all emission interior to Rsqq. 
Right: Quantities estimated in 
the [0.15 - l]/?5oo aperture. 
There is no significant depen- 
dence of 2(r)|bin on temperature 
on a Kendall's r test in either 



In common with most previous work on the subject, we find 
that the slope of the L-T relation of the R EXCESS sample is 
steeper than the prediction from the expectations of self- similar 
collapse models. The steeper slope is found consistently in all 
subsamples, and in all cases the statistical precision of the data 
allow us to rule out the self- similar predictions. We find similar 
results for the L-Yx relation, where the observed slope of ~ 1.0 
is significantly steeper than the self- similar expectation of 0.8, 
and for the L - M500 relation, where the observed slope of 1.8 
is steeper than the expected value of 1.3. These facts imply that 
one or more of the assumptions listed above does not hold for 
the real cluster population. 

The assumption of pure bremsstrahlung is not strictly valid 
since line emission becomes increasingly important as the tem- 
perature decreases, having the eflTect of flattening the relation as 
lower temperature systems are boosted in luminosity. While the 
lower temperature limit of the R EXCESS sample, 2 keV, should 
suflftce to minimise these eflTects, systematic diflTerences in the 
metallicity between objects may serve to change the temperature 
dependence of the X-ray emission from the expected value of 
r^/^. We tested this using the measured temperatures and abun- 
dances of the R EXCESS sample, finding a best fitting power law 
relation of 0.5, in full agreement with the expected dependence. 

The assumption of virial equilibrium leads to the expected 
relation M oc T^^^ between total mass and temperature. A 
topic of vigorous debate in previous years, several recent in- 
vestigations of the X-ray mass-temperature relation have shown 
percent-level agreement in normalisation and that the slop e is not 
greatly different from the self- similar expectation ( Arnaud et al.l 
1200 5: Vikhlinin et al. 2006), although it may be slightly steeper 
dArnaud et al..2005: Sun et al.ii2008.) . 

The question of structural regularity has also received quan- 
titative tests in recent years. For instance, there is now converg- 
ing evidence that the total mass density profiles of galaxy clus- 
ters and groups scale quasi- self- similarly with a mass depen- 
dence that is in good agreement with predictions from numer- 
ical simulations (eg., Pointecouteau et al. 2005; Vikhlinin et al. 
[2006l:lGastaldello et al .120071) . On the mass scales we are consid- 
ering here, the variation of the total mass d ensity concentration 
with mass is in fact consi stent with zero (iPointecouteau et aP 
I2OO5I: lyikhlinin et al.'"2006). This would imply that a variation 
of cluster dark matter structure with mass cannot be responsible 
for the steepening of the L-T relation, at least in the mass range 
covered by the present data. 



However, the baryonic components of clusters are subject to 
somewhat diflTerent physics, and there are indications that the 
large scale ICM density structure is temperature/mass depen- 
dent. The clear correlation of the slope of the gas density profile 
measured in the radial range [0. 3 - 0.817?5oo with te mperature 
seen in the REXCESS analysis o flCroston et al.l (l2008h is one ex- 
ample. At the same timc iCroston et al.' have shown that the tem- 
perature dependence of the relative dispersion of scaled gas den- 
sity profiles has all but disappeared at 0.7 R500, suggesting that 
clusters become increasingly structurally similar at larger radii. 
This is borne out in the present data when we examine the more 
powerfully diagnostic structure factor 2(r)|bin = (Pgas)/(Pgas)^, 
where the average is taken over the radial gas density profile. 
This quantity is identical to that presented in Arnaud & EvrardI 
( 1999) except that we use a fully deconvolved, deprojected gas 
density profiles rather than yS-model fits. 2(r)|bin eflTectively 
probes the variation of the large scale shape of the gas density 
with temperature, that is to say, variation of the gas concentra- 
tion with mass. However, since Pgas(^) is derived from spheri- 
cally symmetric deprojection of the surface brightness profile, 

by construction Pgas(^) is ^(pgas), where (Pgas) is the average 
within each radial shell. Thus we emphasise that 2(r)|bin is only 
a partial estimator of Q(T) which does not probe more subtle 
eflfects such as variations of the gas dumpiness with mass, or 
substructuring at small scale. 

In Figure [7] we plot 2(r)|bin versus system temperature for 
both apertures considered in the present work. The evidence for 
a correlation between 2(r)|bin and system temperature is very 
weak on a Kendall's r test, being significant at only ~ 8 per cent 
for the full aperture and ~ 15 per cent for the [0.15 - l]/?5oo 
aperture. Interestingly, the standard deviation of 2(r)|bin in the 
[0.15 - l]7?5oo aperture, o-Qf^^^ = 0.14, gives an observational 
limit of the variation of cluster structure outside the core regions 
at fixed temperature. That it is only on the order of 15 per cent ar- 
gues for a cluster population of remarkable structural similarity. 
It is thus unlikely that a systematic dependence of cluster struc- 
ture on temperature/mass can be a major cause of steepening of 
the L-T relation of the present sample. However, the eflfect of 
a systematic dependence of gas structure at smaller scale with 
temperature/mass remains an open issue. More detailed assess- 
ment of the gas dumpiness requires combining X-ray data with 
high quality SZ data. 
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Fig. 8. Gas mass fraction vs mass. Left panel: Trend of gas mass fraction versus mass de r ived from X-ray measurements of 41 groups and clusters 
with high quality hydrostatic mass estimates (lArnaud et alj|2007l : IVikhlinin et alj[200^ : ISun et alj|200^ . Grey points are actual measurements; 
black points are mean values in logarithmic mass bins. The solid line is the orthogonal BCES fit to the unbinned data, /gas,500 ^ M^^. Right panel: 
Approximate gas mass fraction versus mass measurements for the REXCESS sample, where the masses have been estimated from the M500 - T 
relation of A rnaud et al.l (1200 5). Black points again show the mean trend for three logarithmic mass bins. The solid line is the same as in the left 
panel. The band illustrates the WMAP 5-Year baryon fraction constraints (iDunklev et al.ll2009h . 



A more likely explanation is if the gas mass-temperature re- 
lation Mgas - T deviates from predictions (or equivalently, if 
the gas mass fraction, /gas = M^^sml^sm varies with mass), 
as was discussed in previous work based on ROSAT data (e.g., 
iNeumann & ArnaudlboQll: iMohr et al.lll999l) Rece nt results on 
the relaxed cluster samp le of Arnaud et al.l (20Q7h and on the 
REXCESS sample itself (ICroston et al."'2008') have shown that 
^ r^, in agreement with previous work, thus implying a 



,500 . 



Steeper dependence than the self- similar prediction of Mgas ^ 
r^/^. Such a dependence of Mgas on T would imply /gas oc Ml/3 
for a self-similar M oc T^^^ relation. 

The left hand panel of Figure[8] shows the /gas,500 - ^500 rela- 
tion for 41 systems ranging from 10^^ to 10^^ M© for which gas 
mass fraction estimat es derived from hydrostatic mass measure- 
ments are avail able (lArnaud et al.1 l2007l: IVikhlinin etaP l2QQ6l: 
ISun et al.ll2QQ8|) . The trend of mean gas fraction with mass is 
evident to the eye. To better illustrate this, we have divided the 
data into five approximately equal logarithmic total mass bins 
and calculated the mean and standard deviation of the gas mass 
fraction measurements in each bin. A BCES orthogonal fit to 
the combined (not binned) data set in log-log space yields the 
relation: 



/i(z)'/2ln/gas,5oo = (-2.37 + 0.03) 

+(0.21 + 0.03) X In (M/2 x 10^^ Mq) 

with (Tinf = 0.12 + 0.03 dispersion. The right hand panel of 
Figure [8] shows the corresponding trend of /gas,500 with clus- 
ter mass for the REXCESS sample, wher e the data valiies hav e 
been estimated using the M-T relation of lArnaud et al.l (l2005i) . 
Dividing the data into three mass bins and averaging the gas 



mass fraction measurements in each bin yields the thick solid 
points, which are in good agreement with the best fitting rela- 
tion to the combined data from hydrostatic estimates discussed 
above. 

The trend in /gas implies a decrease in gas content in 
poorer systems relative to higher mass systems, a fact which 
manifests itself in an increase in ICM entropy and conse- 
quent suppression in luminosity at lower masses. The two 
most likely physical possibilities are a variation with mass 
in the eflftciency of conversion of baryons into stars or 
in situ non-gravitational energy input from e.g., supernova - 
feedback or AGN (e.g. Puchwein, Sijacki, & Springel| l2008l: 
iBower. Mc Carthy, & Benson 2008), or a combination of the 
two. In a forthcoming paper we will use the entropy distributions 
to probe the source(s) and extent of the entropy redistribution. 

Better constraints in the group regime are still required, es- 
pecially in the light o f the i ncrease d scatter ap parently seen there 
( Osmond & Ponmanl[2QQ4l: but see lSun et al][20081 

5. Conclusions 

We have presented a detailed study of the luminosity scaling re- 
lations of REXCESS, a galaxy cluster sample selected by X-ray 
luminosity in such a way as to optimally sample the cluster X- 
ray luminosity function. REXCESS contains objects of diflTerent 
dynamical states with a range of core X-ray properties, allow- 
ing us to investigate the eflTect of the presence of these systems 
on the scaling relations. The homogeneous nature of the sam- 
ple data, which have all been observed with the same satellite 
to approximately the same depth, combined with an analysis ap- 
proach based on extraction of relevant quantities within scaled 
apertures, has been designed to minimise measurement scatter. 
We found the following results: 
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- The slope of the luminosity-temperature, luminosity- and 
luminosity-mass relations are all steeper at greater than 99 
per cent confidence, than expected for self-similar gravita- 
tional collapse scenarios. 

- The dependence of the radially averaged structure factor 
G(^)|bin on temperature, where both quantities have been es- 
timated for the first time within 7?5oo, is not significant either 
when measured from all emission, or from core-excluded 
emission. This suggests that, contrary to previous results, 
structural variation cannot be a significant contributor to the 
steepening of the relations unless there is a very strong tem- 
perature/mass dependence of gas dumpiness. Furthermore, 
the scatter in 2(r)|bin measured in the [0.15 - l]7?5oo aper- 
ture is ~ 15 per cent, illustrating the remarkable structural 
similarity of the present sample in the outer regions. 

- There is strong evidence for a decrease in gas mass con- 
tent in poorer systems relative to higher mass systems, or 
in other words, a dependence of the gas mass fraction on 
total mass. This eff'ect is clearly seen in the M gas - T rela- 
tion of the present sample (ICroston et al.ll2QQ8) and has been 
seen in many previous samples. Using the total mass deter- 
mined from published hydrostatic estimates of a combined 
sample of clusters and groups spanning 10^^-10^^ M©, we 
find that the gas mass fraction depends on total mass such 
that /gas oc The trend in the REXCESS data, when 
total masses are estimated using the M-T scaling relation, is 
similar. This dependence is the dominant cause of the steep- 
ening of the X-ray luminosity scaling relations. 

- For the whole sample, the scatter of the X-ray luminosity- 
temperature relation derived from all emission interior to 
7?5oo is up to 75 per cent depending on the exact fitting 
method; the scatter is less than 40 per cent for the X-ray 
luminosity- 7x relation. The scatter is not strongly compat- 
ible with a Gaussian distribution in either case, the distribu- 
tion being characterised by a tail caused by the presence of 
cool core systems. 

- Cooling core systems, the one-third of the sample with the 
highest central density, describe the high luminosity enve- 
lope and contribute the majority of the variance to the re- 
lations. A fit only to the cool core systems suggests that 
they are off'set from the relation by a simple normalisation 
factor. However, there is nearly two times more logarithmic 
scatter about the X-ray luminosity-temperature relation for 
cool core systems compared to that for non cool core sys- 
tems. This suggests that there are large diff'erences in the core 
structure even for cool core systems as a class. 

- Systems exhibiting morphological substructure tend towards 
the lower luminosity envelope of the relations. Partly this is 
due to the increase in normalisation of the total sample due 
to the presence of cool core systems. The scatter about the X- 
ray luminosity-temperature relation of morphologically dis- 
turbed systems is 65 per cent, compared to a scatter of 48 per 
cent for cool core systems; within the uncertainties, the vari- 
ance is in fact identical. Furthermore, the scatter in morpho- 
logically disturbed systems is identical to that for morpho- 
logically relaxed systems, due partly to the preponderance 
of cool core systems in the relaxed subsample. 

- Simple exclusion of the emission interior to 0.157?5oo re- 
sults in a reduction of scatter in all relations. For the X- 
ray luminosity- temperature relation, the natural logarithmic 
scatter is 30 per cent, a reduction of more than a factor of 
two. Similarly significant reductions are seen in other rela- 
tions. After exclusion of the core, the scatter in luminosity- 
temperature and luminosity- Fx relations is well described 



with a Gaussian distribution at > 85 per cent confidence. A 
reduction in scatter can also be achieved by considering the 
central gas density, Ue^, as a third parameter in the scaling 
relations. The L-T - rie^ relation has a natural logarithmic 
scatter of 47 per cent; the L-Yx- ^e,o relation has a scatter 
of 22 per cent. 

- Using Yx as a mass proxy, a Malmquist bias corrected lumi- 
nosity mass relation for REXCESS is steeper than the raw 
relation due to the under-representation, for a given mass, of 
low luminosity clusters in the sample. 

The behaviour of the observed luminosity scaling relations 
thus appears to be driven principally by a mass dependence of 
the total gas content. Plausible physical explanations for the de- 
pendence are a variation with mass in the efficiency of con- 
version of baryons into stars or in situ heating after accretion. 
Greater understanding of the source of the dependence will re- 
quire deep observations of a similarly representative sample of 
group scale haloes, to measure accurate luminosities and probe 
the underlying physical causes; furthermore, precise calibration 
of the evolution of the scaling relations is needed, ideally with a 
similarly- selected distant cluster sample, to probe the eff'ect over 
time. 
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Appendix A: Cluster image gallery, sorted by cool 
core and (w) classification 

As an aid to visualisation of the sample, in this Appendix we 
present the images of each cluster that were used to calculate 
the centroid shift parameter (w). We remind the reader that (w) 
was evaluated with the central 0.1 7?5oo excised, to avoid biasing 
as a result of the highly peaked surface brightness of cool core 
systems. The images are derived from the three EPIC detectors 
and have been corrected for vignetting; in addition, point sources 
have been removed and replaced by Poisson noise sampled from 
counts in an annulus surrounding the excised source. 

Figures lA.ll and IA.2I show the resulting images. Contours 
increase in steps of ^/2. As in the original REXCESS paper 
(iBohringer ^t al. 2007), the colour table of each panel is scaled 
by a factor of L[0.1 - 2.4]^-^^, derived from the theoretical rela- 
tionships between surface brightness and radius (Sx ^ R), radius 
and mass (R oc M^^^) and luminosity and mass (L oc M"^^^). 

Clusters are divided into cool core and non cool core sub- 
samples and then arranged in order of increasing centroid shift 
parameter (w). Cool core systems generally appear more mor- 
phologically undisturbed; however, the increase in morphologi- 
cal complexity is always evident at higher values of (w). Note 
that two clusters are classified as both cool core and morpho- 
logically disturbed: RXC J1302 -h0230 and RXCJ2319 -7313. 
Figure Em shows that they exhibit both centrally peaked surface 
brightness and morphological complexity, as expected. 



Appendix B: Survey luminosities and associated 
scaling relations 

B.1. Introduction 

Past (e.g., ROSAT) and future (e.g., eROSITA) X-ray survey 
satellites are primarily sensitive in the soft X-ray band. In this 
Appendix we give the soft band luminosities of the REXCESS 
sample, together with the scaling relations between luminosity 
and temperature and Yx. All quantities are calculated as de- 
scribed in Section [T2l 



B.2. Correction for selection bias 

Selection eff'ects can bias the observed scaling relations if not 
all clusters are completely sampled in a well defined test vol- 
ume. The classical problem is the selection bias of an X-ray 
flux-limited sample, where the more luminous clusters are sam- 
pled from a larger volume than the less luminous clusters. 
Thus for any cluster property that is correlated with luminos- 
ity and features a finite scatter, the more luminous clusters will 
always be over-represented with respect to the less luminous 
clusters. This eff'ect has previoi isly been accounted for in flux- 
limited surveys in wo r ks by e -g--llkebe et al.l ( 2002h:IStanek et al.l 
( 2006b : lPacaud et al.1 (l2QQ7l) : I Vikhlinin et al.1 (l2008l) . In the case 
of REXCESS, the survey geometry is more complex, and has 
properties of both a flux- and a volume-limited sample. 

The basis for the calculation of the bias efl'ect is the sur- 
vey selection function, which for our case specifies the size of 
the sampling volume as a function of X-ray luminosity. The se- 
lection function is given by three ingredients (i) the sky area 
covered by REFLEX (~ 4.24 st) and the luminosity intervals 
and redshift shells defined for the REXCESS sample selection, 
(ii) the incompleteness of the sky coverage due to low expo- 
sure regions of the ROSAT All-Sky Survey (RASS), and (iii) the 
flux measurement error which introduces a dispersion between 
the observed and true X-ray luminosity. For the luminosity and 
redshift boundaries we use the definitions given in Table 1 of 
IBohr inger et al.l (l2007h . For those bins containing more clusters 
than the four clusters per bin selected for REXCESS we use an 
appropriate fractional weighting factor. The incompletene ss is 
given by the sensitivity functio n in Table 8 of Bohringer et al.l 
(120041) . illustrated in Fig. 23 of IBohringer et al.1 (l200lk . which 
accounts for the fact that not the complete redshift shell given in 
step (i) is sampled, but there is an incomplete coverage due to 
lack of sensitivity in some areas with low exposure and/or high 
interstellar column density, Nh- An additional factor must be 
used in connection with the sensitivity function since REXCESS 
clusters were selected to have at least 30 counts in the RASS; 
this is taken into account as explained in Bohringer et al. (2001, 
2004). For the first seven bins we used a further restriction in the 
REXCESS selection, Nh < 6 x 10^^ cm"^, so for these bins we 
have recalculated the sky coverage for the restricted sky region 
meeting the Nh criterion. This and some other corrections are 
small in our case (e.g., the part of the sky covered up to the nom- 
inal flux-limit for a minimal detection of 30 photons increases 
from 78% to 79.5% for the low Nh region), and will not signifi- 
cantly aff'ect the final results here, but it is nevertheless important 
to have all these eff'ects under control. 

The resulting selection volume function shown in the left 
hand panel of Figure iBJl is a step function for the 9 luminosity 
bins of REXCESS, with the steps slightly tilted due to correction 
(ii). Folding in the scatter between measured and intrinsic lumi- 
nosity, which is assumed to be ~ 10%, smooths the step function. 
It is monotonically increasing with X-ray luminosity. It implies 
that luminous clusters are more represented than less luminous 
systems within REXCESS. This can then be straightforwardly 
included in a weighting factor for the clusters of diff'erent lumi- 
nosity for a given mass. One approach is to give each cluster 
a weighting factor inverse to the selection volume ratio of the 
measured and nominal luminosity for given mass before fitting a 
scaling relation. Or, as we do here, we can follow a recipe anal- 
ogous to that described by Vikhlinin et al. (2008, Apendix A.2) 
where the mean bias for given mass is calculated via. 



Bias (InLllnLo) = (InL)-lnLo 
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Fig. A.l. Cool core clusters, sorted from top left to bottom right in order of increasing centroid shift parameter, (w). The images are derived from 
the three EPIC detectors and have been corrected for vignetting. Point sources have been removed and replaced by Poisson noi se sampled from 
counts in a surrounding annulus. Contours increase in steps of ^^2. The colour table of each panel is scaled by a factor of L^^^ (see lBohringer et all 
120071 for details). Note that RXC J1302 -h0230 and RXC J2319 -7313 are classified as both cool core and morphologically disturbed. 



£^ A In L p(ln L\ In Lq) VseiCln L) d\nL 
£^ p(\n L\ In Lo) VseiCln L) J In L 

where Lq is the mean (zero scatter) L for given mass, A In L = 
InL - InLo, and /7(lnL|lnLo) characterizes the scatter of L in 
the mass - luminosity relation, assumed to be lognormally dis- 
tributed, and given by the observed scatter in the L - My relation 
(Table[2l). The bias factor is shown as a function of luminosity in 
the right hand panel of Figure lBJl 

A final subtlety is that the luminosities used as the basis of 
the R EXCESS selection, and thus for the bias calculation above, 
are those calculated as in the original REFLEX catalogue. These 
were iteratively calculated in the [0.1-2.4] keV band in the de- 
tection aperture and extrapolated to an assumed radius of 7?2oo, 
and are th us not equivalent to th e luminosities derived in this 
paper (see iBohringer et afll2004l for luminosity calculation de- 
tails; the appropriate REFLEX lumino sities for the R EXCESS 
sample are given in Table 3 of Bohring er et al.l 12007 ). We fit a 
linear relation in log-log space between the present luminosi- 
ties and those from REFLEX, finding Lreflex = 1.15 x ([0.1 - 
2.4] k^dV Lrexcess )^'^^- The bias correction factor for the appro- 
priate luminosity. Bias (lnL| InLo), is then applied to each data 
point in the sample and the relation is refitted. 

The Malmquist bias-corrected bolometric L - My relation is 
shown in the right hand panel of Figure [6l and the correspond- 
ing fitted power law relation is given in Tabled The correction 
steepens the relation somewhat, due to the under-representation, 
at a given mass, of low luminosity clusters in the REXCESS sam- 
ple. 

B.3. Comparison to other results 

The corrected relations for the two survey bands are also given 
in Table IB. 21 The left hand panel of Figure IB. 21 shows the raw 
and corrected [0.1-2.4] keV band relation compared to previ- 
ous determinations from an X-ray hydrostatic analysis assuming 
isothermality (Reiprich & Bohringer 2002) and a stacked weak 
lensing analysis (RykoflF et al 20081 We convert their L - M200 



relations to M500 using a standard NEW model with a concen- 
tration parameter of C500 = 3.2, the average concentration de- 
rived from the total mass p rofiles of the morphologica lly regular 
cluster sample discussed in lPointecouteau et"an (l2005l) . Our cor- 
rected relation has a 25 per cent higher normalisation than that 
of RvkoflT et al. at our fiducial pivot point of 2 x 10^^ M©. The 
iReiprich & Bohringe^ relation has a 6 per cent lower normalisa- 
tion than our corrected relation at the same mass scale. 

The right hand panel of Figure IB. 21 shows the corrected 
[0.5-2] keV band rela tion compared to the results derived by 
IVikhlinin et all (l2008h using the same Malmquist bias correc- 
tion procedure on a sample of clusters observed with Chandra 
(the Chandra Cluster Cosmology Project, CCCP). The agree- 
me nt in normalisa tion is good at low L/M, but at higher L/M 
the lVikhlinin et al.i relation is somewhat below ours (by approx- 
imately 40 per cent at 8 keV, or 8 x 10^"^ M©). 

The REXCESS and CCCP slopes are slightly diflTerent, al- 
though it is important to note that they are in agreement within 
their Icr uncertainties. The bias correction itself does not play 
a part because there is excellent agreement in the magnitude 
of the scatter about the L - M relation from the two samples. 
We use a different M500 - Yx relation to estimate masses, al- 
though in practice the eflTect of this diflTerence will be small since 
our relation is in good agreement with theirs. One partial ex- 
planation could be due to the systematic offset in measurements 
between Chandra and XMM-Newton, in which, at high temper- 
atures, Chandra overestimates the temperature^. Since in both 
cases masses are derived from Yx = Mgas T, this will have the 
effect of boosting the higher mass Chandra points at a given lu- 
minosity, leading to a flatter relation than the one we find here. 
The effect is of order 20 per cent at 8 keV for a mass calculated 
from the M - Yx relation, which alleviates the difference some- 
what. Finally, the samples contain different clusters. The individ- 
ual samples probe slightly different mass ranges as REXCESS 
contains more lower mass systems, while the local CCCP sample 



^ A Comparison of Cluster Temperatures 

Derived from Chandra and XMM-Newton 

http : //cxc . harvard . edu/cal/memos/hrma_memo . pdf | 
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Fig. A.2. Non-cool core clusters, sorted from top left to bottom right in order of increasing centroid shift parameter, (w). 



contains more higher mass systems. In addition, differences in 
the number of cool core systems, and their distribution across the 
mass range, could change the slope. In particular, if R EXCESS 
has more cool core systems at higher mass than CCCP, their 
higher luminosity would make the R EXCESS relation slightly 
steeper. 
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Table B.l. Cluster properties in the soft X-ray band suitable for surveys. Luminosities are given in the [0.1 - 2.4] keV band appropriate for 
ROSAT, and also the [0.5 - 2] keV band. All quantities are calculated assuming CIm = 0.3, Cl^ = 0.7, and ho = 0.7. 
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Columns: (1) Cluster name; (2) z: cluster redshift; (3) T\ \ spectroscopic temperature of the R < R500 region in keV; (4) L [0.1 - 2.4] 1: [0.1 - 2.4] 
keV band luminosity in the R < Rsoo region in units of 10^^ erg s~^ ; (5) L [0.5 - 2]i : [0.5 - 2] keV band luminosity in the R < Rsoo region in units 
of 10"^"^ erg s~^; (6) T2: spectroscopic temperature in the [0.15 - l]i?5oo region in keV; (7) L[0.1 - 2.4]2: [0.1 - 2.4] keV band luminosity in the 
[0.15 - 1] R500 region in units of 10^^ erg s'^ (8) L [0.5 - 2]2: [0.5 - 2] keV band luminosity in the [0.15 - 1] R500 region in units of 10^^ erg s'^ 
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Table B.2. Observed survey band X-ray luminosity scaling relations for the full R EXCESS sample. For each set of observables (L,A), we 
fitted a power law relation of the form h(zTL = C(A/Aoy, with Aq = 5 keV and 2 x 10^"^ M© keV, and w = -1, -9/5 and -7/3 for T, Yx and 
M, respectively. Results are given for the BCES (Y|X) and BCES orthogonal fitting methods (see Section 1241 ). The intrinsic natural logarithmic 
scatter about the best fitting relation in the In-ln plane is given in each case. ^ Since M is derived from Yx, the values of the scatter in the L - M 
relation are identical to those for the L-Yx relation. ^ Relations corrected for Malmquist bias. 
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Li/Ti: luminosity/temperature interior to R500', 

L2/T2: luminosity/temperature in the [0.15-1] Rsoo ape rture; 

My: total mass estimated from the M500 - Yx relation of lArnaud et al.l (l2007h . 
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Fig. B.l. Left panel: Survey volume for the R EXCESS sample. The dashed line is the raw survey volume; the solid line is the volume folded 
with an assumed measurement error of 10 per cent on the luminosity. Right panel: Malmquist bias for the REXCESS sample as a function of 
luminosity. 




Fig. B.2. Left panel: The L[0 . 1 -2.4]- My relation compared with previous determin ations from X-ray hydrostatic analysis assuming isothermality 
(iReiprich & Bohringe3l2002h and stacked weak lensing analysis dRvkoff et al.ll2008h . The points are the bias-corrected REXCESS values. Right 
panel: Malmquist bias corrected L[0.5 - 2] - My relation compared with the results from lVikhlinin et all (120081) . The best fitting relation is given 
in Table IB. 2l All fits have been undertaken with the BCES orthogonal fitting method. 



